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A search for a narrow diphoton mass resonance is presented based on data from 3.0 fb~^ of 
integrated luminosity from pp collisions at i/s = 1.96 TeV collected by the CDF experiment. No 
evidence of a resonance in the diphoton mass spectrum is observed, and upper limits are set on the 
cross section times branching fraction of the resonant state as a function of Higgs boson mass. The 
resulting limits exclude Higgs bosons with masses below 106 GeV /(? at a 95% Bayesian credibility 
level (C.L.) for one fermiophobic benchmark model. 
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The standard model (SM) of particle physics has 
proven to be an extremely robust theory through its ac- 
curate predictions of many experimental results obtained 
over the last few decades. Although the Higgs mecha- 
nism ^ was proposed in the 1960's, the particle it pre- 
dicts, the Higgs boson (h), has yet to be observed in 
nature. 

The SM prediction for the ft, ^ 77 branching fraction is 
extremely small (reaching a maximal value of only about 
0.2% at a Higgs boson mass (to,,) - 120 GeV/c^) '2]; 
however, in "fermiophobic" models, where the coupling 
of the Higgs boson to fermions is suppressed, the dipho- 
ton decay can be greatly enhanced. This phenomenon 
has been shown to arise in a variety of extensions to the 
SM [3,, 4, 5, 6, Jj, and the resulting collider phenomenol- 
ogy has been described 0, 0, [l3| . For this fermiophobic 
case, the diphoton final state dominates at low Higgs bo- 
son masses and is therefore the preferred search channel. 



Fig. [Ija)), and vector boson fusion (abbr. VBF, shown 
inFig.mb)). 

Each of the four experiments [ll|, [ij, llSl, [14 1 at the 





FIG. 1: The dominant production diagrams for the bench- 
mark fermiophobic Higgs boson model: associated production 
with a vector boson (a), and vector boson fusion (b). 



A benchmark fermiophobic model is considered in 
which the Higgs boson does not couple to fermions, 
yet retains its SM couplings to bosons. In this model, 
the fermiophobic Higgs boson production is dominated 
by two processes: associated production (shown in 
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LEP electron-positron collider at CERN place 95% C.L. 
lower limits on the fermiophobic Higgs boson mass (the 
most stringent being 105.5 GeV/c^), while a combina- 
tion of these results obtains a 95% C.L. limit of 109.7 
GeV/c^ [H]. The CDF and DO experiments at the 
Tevatron also searched for a fermiophobic Higgs boson 
[l^ . . Most recently, the DO experiment set limits on 
the production cross section of the fermiophobic Higgs 
boson with 1.1 fb^^ of data, resulting in a 95% C.L. 
lower limit on to,, of 100 GeV/c^ T^. In this Letter, we 
search the diphoton mass distribution from the Collider 
Detector at Fermilab (CDF) for a narrow resonance that 
could reveal the presence of a fermiophobic Higgs boson. 

We use the CDF II detector [13, m] to identify (ID) 
photon candidate events produced in pp collisions at -y/s 
= 1.96 TeV. The innermost detector component is the sil- 
icon vertex tracker [2l| which is surrounded by an open- 
cell drift chamber (COT, [2^). Both sample the trajecto- 
ries of charged particles and determine their momentum 
as they curve in the presence of a 1.4 T axial magnetic 
field. 

Particles that pass through the COT reach the electro- 
magnetic and hadronic calorimeters [i^ 0, [111, which 
are divided into two regions: central {\ri\ < 1.1) and for- 
ward or "plug" (1.1 < \ri\ < 3.6). At the approximate 
electromagnetic shower maximum, the calorimeters con- 
tain fine-grained detectors [2^ that measure the shower 
shape and centroid position in the two dimensions trans- 
verse to the shower development. The calorimeters are 
surrounded by a system of muon chambers . 

Three levels of real-time event selection (trigger) sys- 
tems are used to filter events. The trigger paths used 
here require two clusters of deposited energy in the elec- 
tromagnetic calorimeter. One path requires that both 
clusters have a transverse energy Et > 12 GeV and 
be isolated from other energy clusters in the calorime- 
ter [2^. A second trigger has a cluster transverse en- 
ergy requirement of Et > 18 GeV without the require- 
ment of cluster isolation. By combining these two trigger 
paths, virtually all of the identifiable diphoton events are 
recorded. 

The analysis is divided into two independent subsam- 
ples according to the position of the photons: the first 
requires that both photons be located within the fidu- 
cial region of the central electromagnetic calorimeter 
(jryl < 1.05), and the second requires that one photon be 
located in this region and the other in the plug calorime- 
ter (1.2 < \ri\ < 2.8). The former will be referred to 
as central-central (CC) events, and the latter as central- 
plug (CP) events [29|. The data were recorded between 
February of 2002 and April of 2008, corresponding to an 
integrated luminosity of 3.0 fb~^ for CC and 2.9 fb^^ for 
CP events. 



A series of baseline selection criteria helps to remove 
background events and to ID high-energy photon candi- 
dates for the analysis. Individual photons are required to 
have Et > 15 GeV, while the diphoton pair is required 
to have mass rrij^ > 30 GeV/c^. Photons are required 
to pass CDF standard photon ID requirements including 
the following 28, 3^: transverse shower profiles consis- 
tent with single photon expectation from test beam stud- 
ies [3l| , additional transverse energy in the calorimeter in 
a cone of angular radius R = ^(A0)2 4- (Atj)'^ = 0.4 ^20] 
around the photon candidate be less than 2 GeV, and the 
scalar sum of the px of the tracks in the same cone be 
less than 2 GeV/c. Central photons must also be isolated 
in the shower maximum detector. 

The above selection criteria define an inclusive dipho- 
ton sample. However, the fermiophobic Higgs boson is 
only produced at a non-negligible rate in association with 
a.W 01 Z boson or via the VBF process. In order to im- 
prove sensitivity, the event selection was further extended 
to take advantage of the final state features present in 
these production modes. Associated production domi- 
nates the production process (for mh — 100 GeV/c^ its 
rate is about four times larger than VBF), so the opti- 
mization was carried out on the basis of the associated 
production process alone. A selection based on the fol- 
lowing observables was optimized: diphoton transverse 
momentum {p^'^), transverse momentum of the second 
highest pt jet (p-/^) for hadronic decays of W/Z, and 
missing transverse energy (^t) or transverse momentum 
of an isolated track {p-^''°) for leptonic decays of W/Z. 

For the optimization study, a Bayesian method with a 
flat prior probability was used to estimate the expected 
limits based on signal and background event expectations 
in a 10 GeV/c^ mass window centered at 100 GeV/c^. 
The diphoton background is composed of SM dipho- 
ton events 25 %) and events in which either one or 
both photon candidates are actually quark or gluon jets 
which were misidentified as photons (~ 75 %). Higgs 
boson events with only the diphoton decay mode and 
SM diphoton events were generated using the PYTHIA 
6.2 i52j Monte Carlo event generator and a parametrized 
response of the CDF II detector Issl. Isl . All PYTHIA 
samples were made with CTEQ5L '35'] parton distribu- 
tion functions, where the PYTHIA underlying event model 
is tuned to CDF jet data [s^l. The background compo- 
nent arising from jets misidentified as photons was es- 
timated using photon identification control regions from 
data. The control regions do not overlap with the signal 
region, as the events in the control region are required to 
fail at least one of the standard electromagnetic energy 
fraction or isolation requirements, yet pass a looser set of 
these requirements. 

The optimization shows that a requirement of p^?^ 
> 75 GeV/c is approximately as sensitive as any com- 
bination of the other selection criteria. With this re- 
quirement on Pj?''', roughly 30% of the signal remains 



(slightly varying with m/J while more than 99.5% of the 
background is removed. Although the cut was optimized 
based on associated production, VBF also has a higher 
average p^'' than the background processes and is in- 
cluded in the analysis with the same selection. 

The detector acceptance for signal events is calcu- 
lated using the PYTHIA event generator samples described 
above. Since a pure sample of reconstructed photons is 
not available in the data, corrections to the photon iden- 
tification efficiencies due to imperfections in the detector 
simulation are derived using electrons from Z boson de- 
cays. This is justified since the energy deposition in the 
EM calorimeter by electrons and photons is almost indis- 
tinguishable. The electrons are selected with a slightly 
modified version of the photon ID requirements to allow 
the presence of a high pT track. A correction factor to the 
ID efficiency of the simulation of 0.97 (0.94) is derived for 
central (plug) photons by comparing ID efficiencies from 
the detector simulation with the ID efficiencies measured 
in data. 

The largest systematic uncertainties on the expected 
number of Higgs boson events arise from the luminos- 
ity measurement (6%), varying the parameters control- 
ling the amount of initial and final state radiation from 
the parton shower model of PYTHIA (4%) [33, and the 
PYTHIA modeling of the shape of the p^^ distribution 
for the signal (4%). The latter uncertainty was ob- 
tained by studying the effect on the acceptance from 
the differences in the shape of the pr^^ distribution from 
leading-order, next-to-leading-order, and PYTHIA predic- 
tions [38] . Other systematic uncertainties were also con- 
sidered due to uncertainties in photon ID efficiency, the 
electromagnetic energy scale, and parton distribution 
functions The signal acceptances are included in 

Table H] and they have a relative uncertainty of 8% (9%) 
for CC (CP). 

The decay of a Higgs boson into a diphoton pair ap- 
pears as a very narrow peak in the invariant mass distri- 
bution of these two photons. The diphoton mass resolu- 
tion as determined from simulation is better than 3% for 
the Higgs boson mass region studied here and is limited 
by the energy resolution of the electromagnetic calorime- 
ters. The simulated resolution was compared with data 
using electrons from Z boson decays and is found to 
model the detector well. The diphoton invariant mass 
distribution shown in Fig. [2] could be sensitive to a res- 
onance in the context of the fermiophobic benchmark 
model (see the inserts in Fig. [5] for the signal shape ex- 
pected from simulation). No evidence of such a resonance 
is apparent in the data. As a background prediction 
to be used for setting limits on a diphoton resonance, 
the data is fit to a sum of two exponentials multiplied 
by a fractional degree polynomial, where the degree of 
one term is a parameter of the fit. The fit, excluding 
a 10 GeV/c^ mass window centered at each test m/i, is 
performed for each rah hypothesis (10-GeV/c^ steps from 
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FIG. 2; The invariant mass distribution of central-central (a) 
and central-plug (b) photon pairs after the requirement of 
Pj7^ > 75 GeV/c, shown with the fit to the data for the 
hypothesis of a rrih of 100 GeV/c^. The gap in the fit centered 
at 100 GeV/c^ represents the signal region for this mass point 
that was excluded from the fit. The error bands show the 
statistical uncertainty in the fit. The expected shape of the 
signal from simulation is shown in the inserts. 



70 to 150 GeV/c2). The fit for a nih of 100 GeV/c^ is 
shown in Fig. [2l 

After the background fit for each mass hypothesis has 
been determined, the presence or absence of a Higgs 
boson signal is ascertained on the basis of a binned 
likelihood method incorporating the simulated signal 
shape and the systematic uncertainties. We calculate 
a Bayesian C.L. limit for each mass hypothesis based 
on the combined binned likelihood of the mass distri- 
butions for the CC and CP samples. A posterior density 
in axB{h — > 77) is obtained by multiplying this like- 
lihood by Gaussian prior densities for the background 
normalizations and systematic uncertainties leaving the 
production cross section {axB{h — 77)) with a uniform 
prior density. A 95% C.L. limit is then determined such 
that 95% of the posterior density for axB{h — > 77) falls 
below the limit 4l| . 

The results of the limit calculation are included in Ta- 
ble [T] and displayed graphically in Fig. [3l The SM cross 
sections assumed in the benchmark fermiophobic model 
are used to convert the limits on axB{h — > 77) into lim- 
its on B{h 77). The result excludes the benchmark 
model predictions (at 95% C.L.) for rrih of less than 106 
GeV/c2. 

This Letter presents the results of a search for a nar- 
row resonance in the diphoton mass spectrum using data 
comprising 3.0 fb~^ of integrated luminosity from the 



CDF II detector at the Tevatron. Selected events in- 
clude two photons in the central detector or one photon 
in the central and one in the plug detector. There is no 
evidence of a narrow resonance. Limits are placed on 
the production cross section and the branching fraction 
for the Higgs boson decay into a photon pair and com- 
pared to the predictions of a benchmark fermiophobic 
model resulting in a limit on the Higgs boson mass of ruh 
> 106 GeV/c^ at the 95% C.L. This mass limit is ap- 
proximately as strong as any previous single experiment, 
and the result significantly extends the excluded region 
of B{h 77) for ruh above 110 GeV/c^. 
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FIG. 3: The 95% C.L. upper limit on the branching fraction 
for the fermiophobic Higgs boson decay to diphotons, as a 
function of nih- The shaded regions represent the one and 
two sigma probability of fluctuations of the observed limit 
away from the expected limit based on the distribution of 
possible experimental outcomes. For reference, the 95% C.L. 
limits from LEP are also included. 
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TABLE I: Observed and expected 95% C.L. limits on the production cross section and branching fraction and theory predictions 
for the fermiophobic benchmark Higgs boson model. Total signal acceptance is also included. Note that the CP channel is not 
included for the 90 GeV/c^ point 29]. 
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